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Abstract 

The QCD analysis of the xF 3 structure function measured in deep-inelastic 
scattering of neutrinos and antineutrinos on an iron target at the Fermilab 
Tevatron is done in 1-, 2- and 3-loop order of QCD. The x dependence 
of the higher-twist contribution is evaluated. The experimental value of 
higher-twist corrections to the Gross-Llewellyn Smith sum rule is discussed. 



At present, the precise measurements of structure functions (SF) and detailed the- 
oretical calculations of QCD predictions for scaling violations ( up to 3-loop order for 
xF 3 (x, Q 2 ) and F 2 (x, Q 2 )) provide an important means of accurate comparison of QCD 
with experiment. The importance of higher-twist (HT) contribution to SF was pointed 
from the very beginning of QCD comparison with experimental data [J] on SF. De- 
spite a fast progress in theoretical QCD calculations of power corrections to nonsinglet 
SF and sum rules ||, [J ( for reviews and references see [f|]), the shape of HT (order 
1/Q 2 ) contributions is measured only for F 2 SF |5| and is still only estimated for xF 3 
0. In the present note, the x dependence of HT contribution is phenomenologycally 
determined in the framework of QCD analysis of the experimental data of the CCFR 
collaboration obtained at Fermilab Tevatron |7j for the xF% structure functions of the 
deep-inelastic scattering of neutrinos and antineutrinos on an Iron target by means of 
the Jacobi polynomial expansion method in the 1-, 2- and 3-loop order of QCD. 

The details of this method are described in ||-[]13||. The Q 2 - evolution of the 
moments M 3 QCD (N,Q 2 ) is given by perturbative QCD [fH|, |T5|. 
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Here a s (Q 2 ) is the constant of strong interaction, r j ( ^ >NS are nonsinglet leading order 
anomalous dimensions. The factor (Qq, Q 2 ) contains all next- and next-to-next- 
to-leading order QCD corrections and is constructed in accordance with [pl| based on 
theoretical results of (|TjJ. 

Having at hand the moments (|l|) and following the method ||, 10 , we can write the 
structure function xF 3 in the form: 
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where Q^(x) is a set of Jacobi polynomials and c"(a,/3) are coefficients of the series 
of 0"'^(x) in powers of x: 

n 

^ n {x) = Y,cf\P)x>. (3) 

j=0 

The unknown coefficients M 3 (N, Qq) in ([!]) could be parametrised as Mellin moments 
of some function: 



M 3 QCD (N,Q 2 ) = [ dxx N - 2 Ax b (l -x) c (l + -fx), N = 2,3, ... 
Jo 



(4) 



For N„ 



8 the accuracy better than 10 3 is achieved in a wide region of 



parameters a and (3 [EX]. In particular, we use a = 0.7 and (3 = 3.0 



1 Announced by CCFR collaboration reevaluation of the structure functions could change the results 
of the QCD analysis. 
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Using Mellin moments (0),(§), expression (|2[) for SF and taking target-mass correc- 
tions (TMC) into account, we have reconstructed xF^ CD \x, Q 2 ). Five free parameters: 
A, b, c, 7 and QCD parameter are to be determine from comparison with experi- 
mental data. 

To extract the HT, contribution we parameterize the nonsinglet SF as follows: 

xF 3 (x, Q 2 ) = xFi QCD (x, Q 2 ) + h(x)/Q 2 , (5) 

where the Q 2 dependence of the first term in the r.h.s is determined by perturbative 
QCD. Constants h(x,j) (one per x-bin) parameterize the HT x dependence. In accor- 
dance with the x-bin structure of the CCFR data we put Xj = 0.015, 0.045, 0.080, 
0.125, 0.175, 0.225, 0.275, 0.350, 0.450, 0.550, 0.650 for i = 1, 2.. .11. The values of 
constants h(xt) as well as parameters A, b, c,7 and scale parameter A are determined by 
fitting the set of the CCFR data at 90 experimental points of xF 3 in a wide kinematical 
region: 1.3 GeV 2 < Q 2 < 501 GeV 2 and 0.015 < x < 0.65 and Q 2 = 10 GeV 2 . 
We have put the number of flavours to equal 4. The TMC are taken into account to the 
order of o(M^ ucl /Q A ) . The nuclear effect of the relativistic Fermi motion is estimated 
^from below by the ratio R^ N = F^/F^ []TH| obtained in the covariant approach in 
light-cone variables fi~7|l . 

Results of the fit are presented in Table 1 and Figures 1-3. The theoretical prediction 
for h(x) from is presented at Figure 3. 



Several comments are in order: 

• A decrease of ■^ 2 ( NNL °) in comparison with ^ 2 ( 7VLO ) and ^ 2 ( L °); Xd!f NL < Xdj < 

2 LO 

Xjf demonstrates that 3-loop effects are important for the kinematical region 
under consideration. For all orders of QCD the \ 2 P er degree of freedom is smaller 



than in where the fit was done without HT contribution. 



The obtained value of the A is smaller in comparison with results of the previous 
analysis of CCFR data 0, with the cut off Q 2 > 10 GeV 2 A^ L) = 184 ± 
31 MeV but exhibits relatively large statistical errors. Results of the NNLO 
fit gives the constant of strong interaction ag NLO (M^) = O.lOAt Q^lisyst.) in 
agreement within the errors with usual DIS results Jl9| and with the predictions 
of CCFR-NuTeV Collaboration based on the test of the Gross-Llewellyn Smith 
(GLS) sum rule. 

The shape of h(x) demonstrates for LO, NLO and NNLO fit a very small value at 
0.015 < x < 0.045, a negative value at 0.1 < x < 0.045 (with a minimum located at 
about x = 0.2) and increase from a negative to a positive value at 0.2 < x < 0.65. 
This behavior is in qualitative agreement with theoretical predictions of and 
reproduces appropriately the predicted zero of h(x): x theor ~ 0.67 while in our 
NNLO analysis x NNL ~ 0.40 . A separate fit with cuts off Q 2 > 5 GeV 2 and 
Q 2 > 10 GeV 2 shows the stability of shape of h(x) and increase of errors. 
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The absolute value of h(x) slightly decreases from LO to NNLO fit. It may be 
indicates a special role of higher order perturbative QCD corrections reveals by 
renormalon technique p5 |: at higher order xF% D in (|5]) describes effectively the 
power corrections. 

Definite theoretical predictions are presented for the first moment of h(x) which 
contributes to the GLS sum rule h\ = Jq ^-dx . A general structure of this 
contribution is known ^from the results of Ref.[^] The corresponding numerical 
calculations of this term was made in Ref. [[SJ hi = —0.29 ± 0.140 and more 
recently in Ref. |3| hi = —0.47 ± 0.04, using the same three-point function 
QCD-sum-rules technique. One can estimate hi based on the results of Table 1. : 
h{° = 0.12 ± 0.53 , h? LO = 0.14 ± 0.53 and h? NLO = 0.13 ± 0.45. Taking into 
account the errors the values of h± , h^ LO and h^ NLO could be compared with 



the prediction of [23[ and the recent result of p5| for GLS sum rule 
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It should be noted that the fit without the nuclear effect R 



NLO,R=l 



D/N 

F 
1 



1 provides 

h LO,R=i = o.ii ± 0.51 , /if LOA=1 - 0.12 ± 0.40 and hi NLO ' R=l = 0.12 ± 0.48 in 
a good agreement with previous results. The large contribution of small x region 
to hi needs the shadowing correction taking into account for more detail analysis 



In conclusion it should be stressed, that for precise determination of the HT contribu- 
tion to SF the role of nuclear effect should be clarified and a more realistic approximation 



for R 



Fe/N 



Ft 



IF 
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is needed. A possible interplay of the nuclear effect and TMC 
was considered in [ p7|] . We also did not take into account the threshold effects on Q 2 
evolution of SF due to heavy quarks [28[] which is necessary owing to a wide kinematical 
region of data under consideration. 
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Figure captions. 



Fig.l. Higher-twist contributions from LO fit and the theoretical prediction for h(x) 
from [Fl. 



Fig. 2. Higher-twist contributions from NLO fit. 
Fig.3. Higher-twist contributions from NNLO fit. 
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Table I. Results of 1-, 2- and 3- order QCD fit (with TMC) of the CCFR xF 3 SF 
data for / = 4, Q 2 > 1.3GeV 2 with the corresponding statistical errors and values of 
h(x) at different values of x. N MAX = 10 for 1- and 2- oder and N MAX = 7 for 3- order 
fit. 



7 



Table I. 





LO 


NLO 


NNLO 




65.1/74 


62.9/74 


60.9/74 


A 


6.69 ± 0.87 


6.04 ± 0.51 


5.56 ± 0.18 


b 


0.772 ± 0.040 


0.745 ± 0.026 


0.719 ± 0.011 


c 


4.04 ± 0.16 


3.97 ± 0.14 


3.91 ± 0.12 


7 


0.424 ± 0.53 


0.603 ± 0.317 


0.707 ± 0.055 




76 ± 62 


132 ± 80 


134 ± 57 


[MeV] 










h( Xi ) [GeV 2 ] 


0.015 


0.012 ± 0.034 


0.018 ± 0.036 


-0.015 ± 0.022 


0.045 


-0.008 ± 0.049 


0.037 ± 0.063 


0.043 ± 0.054 


0.080 


-0.199 ± 0.061 


-0.107 ± 0.079 


-0.067 ± 0.077 


0.125 


-0.318 ± 0.084 


-0.203 ± 0.083 


-0.144 ± 0.086 


0.175 


-0.175 ± 0.133 


-0.073 ± 0.114 


-0.005 ± 0.106 


0.225 


-0.242 ± 0.186 


-0.176 ± 0.159 


-0.113 ± 0.133 


0.275 


-0.217 ± 0.241 


-0.202 ± 0.210 


-0.162 ± 0.168 


0.350 


0.095 ± 0.294 


0.023 ± 0.253 


0.011 ± 0.185 


0.450 


0.129 ± 0.302 


-0.010 ± 0.280 


-0.051 ± 0.207 


0.550 


0.283 ± 0.235 


0.150 ± 0.249 


0.086 ± 0.205 


0.650 


0.510 ± 0.155 


0.412 ± 0.180 


0.349 ± 0.159 
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